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Abstract 

An exact reduced-density-operator for the output quantum states in time- 
convolutionless form was derived by solving the quantum Liouville equation 
which governs the dynamics of a noisy quantum channel by using a projection 
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and the mutual interaction between the system and the reservoir are given. 
Then, we apply the formulation to model a two-bit quantum gate composed 
of coupled spin systems in which the Heisenberg coupling is controlled by the 
tunneling barrier between neighboring quantum dots. Gate Characteristics 
including the entropy, fidelity, and purity are calculated numerically for both 
mixed and entangled initial states. 



Typeset using REVTgX 



*Also with Department of Electrical Engineering, University of Seoul, Seoul 130-743, Korea 
^Electronic address; dahn@uoscc.uos.ac.kr 

^Permanent address: Department of Electronics Engineering, Korea University, 5-1 Anam, 
Sungbook-ku, Seoul 136-701, Korea 

1 



I. INTRODUCTION 



There has been a considerable interest in the quantum theory of information and com- 
putation for the past several years |j]-|l6j]. Especially, quantum-mechanical properties of 
coding [HU, noisy-channels including error- correcting codes and channel fidelity ||, 



and computation [p!0| --[T5|] have been studied in detail. It was shown [|T3| , |l^| that any quantum 
computation procedure can be decomposed into operations on single-bit gates and a two-bit 
gate which involves an entanglement operation on two quantum bits or qubits. Presence of 
decoherence and imperfections cause the operations of these quantum gates away from the 
ideal ones and as a result one can regard these gates as a part of noisy quantum channels. 
Detailed analysis of these channels are necessary for the complete understanding of general 
quantum information process. Mathematically, the dynamics of quantum channels or gen- 
eralized quantum gates involves the transformation of input quantum states represented by 
a density operator p into an output states p' [fTE], i.e., 



p-^p' = £[p], (1) 

where we assume S is a linear mapping but is not necessarily a unitary transformation if one 
considers an open system interacting with the reservoir such as noisy quantum channels. 
A model of a noisy quantum channel would involve several Hamiltonians for the system 
representing qubits, reservoir and the mutual interaction between the system and the reser- 
voir that causes the decoherence or noise. The density operator is then governed by the 



quantum Liouville equation [17| which is an integro-differential equations and in general, it 
is nontrivial to obtain the solution of the form given by Eq. (|l|). Rather, one is expected 
to get the solution for the density operator for the output states in Volterra type integral 
equation: 

p(t) = A(t, 0)p(0) + J drB(t, r)p(r) (2) 

where A is a propagator and B is a memory kernel. In general, it is very difficult to solve for 
the memory kernels of the time-convolution form equation (§) self-consistently and almost 



always, one must be content with the narrowing limit or the fast modulation limit p0| 



Some time ago, the time-convolutionless equations of motion in the Heisenberg picture 



was suggested by Tokuyama and Mori [0 to overcome above mentioned difficulties for 
problems in nonequilibrium statistical mechanics. These formulations were then developed 
in the Schrodinger picture by using the projection operator technique JT§j-F2T|| . One of the 
authors applied the time-convolutionless formulation to the model of quantum devices for 
detailed numerical study [p^ -p5|. It was shown that the time-convolutionless formulation 
can also incorporate both non-Markovian relaxation and renormalization of the memory 
effects. 



Recently, Loss and DiVincenzo [yj has made a comprehensive study of the two-bit 
quantum gate taking into account the effect of decoherence on the gate operation using the 
reduced density operator in the time-convolution formulation. Their results indicate that 
the detailed analysis of the decoherence process is important for the reliable operation of 
quantum gates utilizing controlled, nonequilibrium time evolution of solid-state spin systems. 

In order to make the reduce-density operator for the output quantum states of the form 
given by the equation ([I]), several approximations including the Born approximation were 
made in their theory. In our opinion, it would be more convenient if there is a way to get 
exact solution for the output density-operator in time-convolutionless form given by ([!]). 

In this paper, we first derive the exact solution for the reduced-density-operator of the 
output quantum states in time-convolutionless form by solving the quantum Liouville equa- 
tion for a quantum channel using the projection operator method. The formalism we develop 
would be general enough to model a realistic quantum channel or a quantum gate. Secondly, 
we apply the theory to model a two-bit quantum gate composed of coupled spin systems in 
which the Heisenberg coupling is controlled by the tunneling barrier between neighboring 
single electron quantum dots. 
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II. TIME-CONVOLUTIONLESS REDUCED-DENSITY-OPERATOR THEORY 
OF A QUANTUM SYSTEM INTERACTING WITH A RESERVOIR 



In this section, we study the quantum Liouville equation for a quantum system which 
corresponds to a quantum channel or a generalized quantum gate to derive an equation and 
to solve for a reduced-density-operator of a system coupled to a reservoir. An interaction 
between the system and the reservoir leads to decoherence. The Hamiltonian of the total 
system is assumed to be 

H T {t) = H s (t) +H B + H int , (3) 

where Hs(t) is the Hamiltonian of the system representing a quantum gate (or channel), H B 
the reservoir and H int the Hamiltonian for the interaction of the system with its reservoir. 
The evolution of the system might include a coding, transmission and decoding process. The 
equation of motion for the density operator pr{t) of the total system is given by a quantum 
Liouville equation 

j t p T {t) = -i[H T ,p T ] 

= -iL T p T , (4) 



where 



L T (t) = L s (t) + L B + L 



hit 



is the Liouville superoperator in one-to-one correspondence with the Hamiltonian. In this 
work, we use a unit where h = 1. In order to derive an equation and to solve for a system 
alone, it is convenient to use the projection operators p6| , |27| which decompose the total 



system by eliminating the degrees of freedom for the reservoir. We define time-independent 



projection operators P and Q as [19 



PX = p B ti B (X), Q = l-P, (5) 
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for any dynamical variable X. Here tre indicates a partial trace over the quantum reservoir. 
Projection operators satisfy the operator identity P 2 = P_,Q 2 = Q and PQ = QP = 0. The 
information of the system is then contained in the reduced density operator p(t) which is 
defined by 

p(t) = tr BPr(t) 

= ti B Pp T {t). (6) 

In order to derive a time-convolutionless equation, we first multiply Eq. (f|) by P and Q to 
obtain coupled equation for P_pr(t) and Qpr(t): 

j t Pp T (t) = -iPp T Pp T (t) + iPL T (t)Qp T (t), (7) 

j t QpAt) = -iQprQpAt) + iQL T (t)Pp T (t). (8) 

We assume that the channel was turned on at t = and the input state prepared at t — 0, 
p(t = 0) was isolated with the reservoir at t — 0, i.e., Qpt(0) = . 
The formal solution of (§) is given by |E| 

Qp T (t) = -i f drH(t, t)QL t (t)P Pt (t), (9) 
— Jo — 

where the projected propagator Hit, r) of the total system is given by 

H(t,r) =Texpl-ij dsQL T (s)Q^ . (10) 

Here T denotes the time-ordering operator. Because Eq. (||) is in time-convolution form, we 
transform the memory kernel in (^) into time-convolutionless form |22[ by substituting the 
formal solution of (f|) 

p T (r)=G(t 1 r)p T (t) (11) 
into Eq. (§). The anti-time evolution operator G(t, r) of the total system is defined by 

G_{t, r) = T c exp Si J* dsL T {s) J , 



where T c is the anti-time-ordering operator. From Eq. (|9|) and (pTTj) , we obtain 

^ = {^-1}^) (12) 

where 

e-\t) = g {t) 

= 1 + % f drH(t, t)QL t (t)P G(tr) (13) 

J 

By substituting Eq. fll2f) into (0), we obtain the time-convolutionless equation of motion for 
Pprit) as 

j t Ep T (t) = -iPL T (t)Pp T (t) - iPL T {t){6{t) - l}Pp T (t) (14) 



It can be shown that the formal solution of ([14]) is given by 

RpT = U(t, 0)Pp T (0) - % [* dsU{t, s)PL T (s){6(s) - l}Pp T {s), (15) 

Jo 

where the projected propagator U_(t, r) of the system is defined by 

U(t, t) = Texp j-i J* dsPL T (s)P^ . (16) 
To transform Eq. ( |15|) into time-convolutionless form once again, we substitute 

p T (s)=G(t,s)p T (t) (17) 

into dl5|) to obtain: 

Pp T (t) = 0)Pp r (0) - i /* dsC/(t, s)PL T (s){9(s) - 1}P G(t, s)Pr(t) 

Jo 

= U(t, 0)Pp r (0) - z / ds£/;(t, s)PL T (s){^(s) - 1}P G(t, s)Pp T (t) 
Jo 

-i ( dsUCt, s)PL T (s){6(s) - 1}P G(t, s)Qp T (t) 
Jo — 

= U(t, 0)Pp T (0) - i I dsUit, s)PL T (s){6(s) - 1}P Git, s)9(t)Pp T (t). (18) 
Jo 

By the way, 



Ppr(t) = p B tr B (f>r(t)) 
= PBp(t), 



and 



Then 



PL T {t)P = P{L s {t) + L B + L mt )P 
= PL s {t)P 
= L s (t)P. 

LL{t, 0)Pp r (0) = Texp j-i jf* dsPL T (s)P j Pp T (0) 
= Pexp | -i jf* rfsL 5 (s)p| Pp T (0) 
= Us(t,0)Pp T (0) 
= Us(t,0)p B p(t). 

Here 0) denotes the propagator of the system. Likewise, 

U(t, s)PL T (s){6(s) - 1}P G(t, s)6(t)Pp T (t) 
= LL s (t, s)p B ti B [L T (s){9(s) - l}p B tT B {G(t, s)6(t)p B }] pit) 
= U s (t,s)p B tT B [L T (s){6(s) - l}p B ]tr B [G(t, s)6(t)p B ]p(t). 

Substituting (0) and (gD into ©, we obtain 

p(t)=LLs(t,0)p(0) 

-i /"* dsU s (t,s)tT B \L T (s){9(s) - l}p B \tT B \G{t,s)e(t)p B \p 

JO 1 J L J 



or 



p(t) = £(t)p(0) 

= W-\t)U s (t,0)p(0), 



with 



Wit) = l + i dsU s (t,s)tT B L T (s){6(s) - l}p B tr B G(t,s)6(t)p B 



Here, we define 



l + i / ds[/ g (t,s)tr B [L illt E(s){l-E(s)}- 1 pB 



xtr 



s(t) = i-r 1 (t), 



(25) 



(26) 



itLg 



U s (t), 



(27) 



and 



R(t, r) = T c exp |i £ dsU^\s)L iat Uv} , 



(28) 



where U_o{t) is the evolution operator of the system with the reservoir and R(t,r) is the 
evolution operator of the total system in the interacting picture. In (p5[), we use the 
identities PL T (s)Q = PL int Q and H(t, r)Q = Q H(t, r). 
Detailed expression for £(t) becomes 

£(t) = 1 - r : (t) 

= -i / drHCt, t)QL t (t)P G(t, r) 
Jo — 

= -i f drH{t, T)QL mt (r)P G(t, r) 
Jo — 

= -i f dTU^(t)S{t,r)U Q l QL mt PU^(r)R{t,r)U^\t), 
Jo — 



(29) 



with 



S(t, t) = Texp |-i J* ds QU- \s)L int U (s)Q 



(30) 

where S(t, r) is the projected propagator |21|] of the total system in the interaction picture. 
It is now obvious from fl2l|) and (|25"D, the exact solution p(t) for the output quantum state 
is in time-convolutionless form given by Eq. ([I]) which is employed in the description of 



quantum information processing and computation [JL6] 



We now consider the case when the system is interacting weakly with the reservoir 
and expand up to the second order in powers of the interaction Hamiltonian H int . The 
renormalization of the unperturbed energy of the system and the first order of the interaction 



H int gives [19-21 



PL mt P = 0. 



(31) 



Then in the lowest order Born approximation which is valid up to the order (ifi nt ) 2 , we 
obtain 



W {2 \t) = l+i dsUs(t,s)tr B L jnb E^(s)p B tr B Uo{s)Uo\t)p B 
Jo 1 J L 

= l + i / dsU s (t,s)ti B \L int i: {1) (s)p B ]Us 1 (t J s), 



(32) 



or 



W {2 \t)] 1 = l-i [ t dsU s (t,s)tr B \L int rt 1 \s)p B ]u s \t,s 
o 



(33) 



and 



£(2) 



l-i / dsU s (t,s)tr B \L int Z (1 \s)p B ]Us\t,s) 
Jo 1 J 



u s (t,o) 



(34) 



Here 



= / drlJ^(s)U^\r)QL int PU J0 (T)U \s) 
Jo — 

= -i dTUj (s,T)L int U^ 1 (s,T). 
Jo 

The time-convolutionless form of the output reduced-density-operator 



(35) 



(36) 

together with (|32])-(p5]) can be used in any time scale and is valid up to the second order in 
powers in the interaction between the system and the reservoir. 

In the next section, reduced-density-operator for the output quantum state is used to 
study the two-bit quantum gate utilizing coupled spin system in nonequilibrium situation. 



III. DECOHERENCE OF TWO-BIT QUANTUM GATE 



We consider a two-bit quantum gate based on nonequilibrium dynamics of the spin of 
excess electrons in quantum dots [[OJ. In this system, the gate operation is controlled 
by an electrical tunneling between two quantum dots. Projecting out the spatial parts of 
wavefunctions of electrons, we model the system by the Hubbard Hamiltonian |H| ; 



H s {t) = J(t)Sx • S 2 (37) 

where J(t) is time-dependent Heisenberg coupling which involves the energy difference be- 
tween the spin singlet and triplet states. If we turn on J(t) for J dtJ(t) = Jqt s = ir, the 
unitary operator associated with the Hamiltonian fl37f) gives the swap operation up to overall 
phase difference; if labels the spin states of two electrons in the S z basis with i,j =f, j, 
then swap operation £/ swap on two registers gives U swap \i,j) = 

In reality, quantum-dot system of our interest is not a closed system, so we have to take 
into account of the decoherence effects due to the interaction with the environment which 
is coupled with the system. For the action of the environment during the gate operation, 



we use a Calderia-Leggett-type model [T5| where a set of harmonic oscillators are coupled 
linearly to the system spins by 

H- mt = A(Si ■ b[ + S 2 ■ b 2 ) (38) 

Here, b{ = J2a9a{ci'ai + a? a i) ls a fluctuating quantum field whose unperturbed motion is 
governed by the harmonic-oscillator Hamiltonian, 

HB{t) = ^u a ala a (39) 

a 

where a\ (a a ) are bosonic creation (annihilation) operator and u a are the corresponding 
frequencies with spectral distribution function A{uj) = n J2 a Qa^i^ — ^a)- 

For a coupled spin system, the evolution operator £^ given by Eq. ( p4|) can be written 
down explicitly in terms of spin operators. Substituting (]3"T|)-(|3Tj|) into definitions for ?7 and 
Li n t, the integrand of Eq. (0) can be written as, 
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tr B [L int X (1) (s)p B \ U s \t, s)Us(t, O)p(O) 

= -i [ drtT B \L- mt U (s, r)L iDti U^ 1 (s, r)p B } U s (s, O)p(O) 

= - Z A 2 [ dr{[Sl S l k (r - s)(U s (s, 0)p(0))]tr B {^L(T - s)p B } 
ijkl J0 

+ [(U s (s, 0)p(0))S l k (r - s),S!}tT B {b l k (r - s)b\p B }} 
= -iJ2f dr{\S!, Sl{r - s)(Us(s, 0)p(0))]{r(r -a)- *A(r - «)} 
+ [(U s (s, 0) P (0))S?(t - s), S!]{T(r - s) + zA(r -«)}}, 
where the trace over the heat bath is done for the harmonic oscillator eigenstates, 



(40) 



Tr B {b l k (t)biPB} = 8 ik 5 fl - I " A(u>) (e-*-* + ^i^. > ^. 



7T JO 



(41) 



and we define T(t) and A(t) as 



T(t) +iA(t) = \ 2 Tr B {bi(t)blp B }. 

Then, Eq. fl34|) leads to 

^ = Jk(t, 0) fl - f ds f dr J2{[S{(s), Si(r)p(0)]{r(r - s) - *A(r - s)} 
+ [p(0)Sf (r), Sf (a)]{r(r - a) + *A(r - s)}} . 



(42) 



(43) 



Now we evaluate the density operator in basis representation; pit) = J2 a p Pap(t)e a p, e a p 
is the basis for the density operators, and in this work we choose e a p as the multiplet states, 
i.e. e Q p = \at)(p\ with a, (3 = 1,2,3,4; |1) = | TT), |2) = (| TI) + I IT»/V2, |3> = | ||), and 
I 4 ) = (I TI) - I IT>)/v / 2. By defining the inner product like (e a p, e 7<5 ) = tr[e^e 7<5 ] = S a p5 lS , 
Pap{t) is obtained as; 



p a p = (e a p,p(t)) 



(e a p,£®p(0)) = Y,(e a p,£ {2) e, 5 )p(0), 5 = £ £^(0)^ 

■yS jS 



(44) 



where p(0) 7 5 expansion coefficients of the initial density operator. Without the interaction 



with environment, i.e. the absence of the second term in Eq. (0), S^phs * s reduced to 
Us(t) a pi s and evaluated on the multiplet basis as 
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(e a p, U s (t)e rS ) = 5 a p5 lS e-^ E ^\ (45) 

where -^1,2,3 = Jo/4 and E4 = —3 Jo/4 are the triplet and singlet energy eigenvalues. Here, 
t has its value t {r s ) if t is less(larger) than r s . Then, U_s(t) becomes the swap operator, 
U s (t)=e- i " / 4 f/ swap ift = r s . 

In order to evaluate 8^ 2 \ we first calculate the following matrix elements: 



(e Q/3 ,E^(^) ) ^( 7 -)^])=E{(«l^(^)^Ml7)(5|/3) 

- (a|^(r)| 7 )(5|^( S )|/3)} 

= ^ E M aKKy e^ +i ^- - M ajSp e if ^ +i ^ (46) 



and 



(e^,EMW.'S?(*)])=E{<*|5?W'S?WI^<«l7> 

- (a|5?( 8 )|7>(tf|5?(r)|/3)} 

= 5 Q7 E M 5 «^e^+^ - M^e*™^^- (47) 
where M a/37(5 = J2ij{a\Si \i3){j\Sl \5), uj a(3 = E a Ep, and 



r s 



t(s) if r(s) < r s 
t s otherwise. 



Then, the matrix element of the evolution operator, S^p^s * s obtained by substituting ( [45]) 
(0) into (H), 



p(2) _ -itui al 



^a 7 ^/3<5 — <5/M E J ^aKK 7 PKK|7a(^) ~ E ^<$«k#P5k|7/?(*) 

K K 

(*)} (48) 



with the time-dependent term p a /3\js{t) defined by 

PaBhsit) = f ' dse~ rsu ^ /"dre i ™^{r(r - s) - iA{t - s)}. (49) 
Jo Jo 

For numerical calculations, it is more convenient to split the time integrals of the matrix 
Pa/3\-ys(t) into three parts; 

12 





dTe iTUJai { 


JO 








f dse' iTsU ^ 5 


[ Ta dre iTUJ 




lo 






f dse- iTsUJ ? s 


f dre iT ^ 




IT a 



= r dse^'r*"*^ f dre iTUJ -< a {T(T) + iA(r)} 
Jo Jo 

+e iT ^ f dse is ^ f dTe iT ^ a {T(r) + zA(r)} 

+e ir s (^+^ 7 ) I dg f s Ta d T {T{T) + iA{t)}. (50) 
Jt s Jo 

In order to investigate the dynamics of the density operator in non-equilibrium situation, we 
calculate Eqs. numerically, assuming an Ohmic damping for spectral distribution 

function A(u) = rju with a cutoff frequency uj c [BU . 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

We now study dynamics of the density operator for various initial states. First, we 
calculate the evolution of the spin states during the swap gate operation and compare 
our results with those obtained by Loss and DiVincenzo |T5] . The initial spin state is 
chosen to be the spin-up for the second electron while the first electron is unpolarized; 
P(0) = (I TT)(TT I + I IT) (IT D/2- In the multiplet basis, the initial state is expanded as; 

P(0) = §I1><1| + \\2)(2\ - \\2){A\ - i|4)(2| + i|4)(4|. (51) 

Fig. l-(a) shows the spin polarization calculated using parameters X 2 r] = 1.8 x 10~ 5 , 
k B T = 300 K, u c = 400 K, and J Q = 1 K(solid lines). For the interval, < t < t s the 
spin polarization of the first electron s = 2 {Si) = 2tr[p(i)Sl] changes to nearly a unity 
whereas the spin state of the second electron becomes zero(dashed line), demonstrating the 
feasibility of the swap operation. However, due to the decoherence, we find that a perfect 
swap operation cannot be achievable. In addition, the perturbing fields cause the monotonic 
decreases of the spin polarization with the elapse time after completion of swap operation. 
This means that spin states are becoming thermalized owing to the interaction with the 
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environment, which shows the decoherence of the states. The decoherence would be a 
fundamental problem in making a reliable quantum logic gate, which puts severe restriction 
on building the realistic quantum computer. However, there are several quantum error- 
correction techniques which can compensate imperfections introduced by the decoherence 
during and after the gate operation Comparing with the result obtained in the previous 

work |15| (dotted line), we find that both calculations yield similar results for t > t s except 
for the value at t = r a . We think that the discrepancy at t = t s is resulted from somewhat 
simplified evaluation of the evolution operator in the reference |H| when the swap operation 
occurs. 

In Fig. l-(b) and (c), we plot the gate fidelity T and gate purity V which characterize 



the intrinsic properties of the gate, and are defined as [^] 



a a,0 



1 



V = tr[p(t)P = 24 E 



|c(2) \2,X^fcP) c(2)* , |c(2) |2 
l°a/3|77l ' \°Q/3|77°a/3|6<5 ' \^aP\f5 

5 



(52) 
(53) 



where the overbar means an average over all possible initial state \tpo) and Us(t) is an ideal 
gate operation which was turned on during the time interval, < t < r s . The last equalities 
in ( |52|) and (|53|) were derived under the condition of both trace and hermiticity of £^ 
being preserved within our approximation scheme. For an ideal quantum gate, the gate 
fidelity T and the gate purity V must be equal to one during the gate operation because 
in that case the evolution operator is unitary. Our calculation shows that both T and V 
are found to decrease almostly linearly as time elapses, which indicates clearly the presence 
of decoherence effect. As the case of the spin polarization, the decreasing rates for T and 



V are close to those obtained in Ref. fTS] , however its value at t = t s are different and our 
results show more severe decoherence of the spin state for the same parameters. 

Another interesting property of the two-bit gates is the von Neumann entropy A defined 
as A = — tr[p(i) log 2 p(t)] of a quantum state. In Fig. l-(d), the calculated von Neumann 
entropy of the spin system is plotted. For the initial density operator of Eq. (|5l|), its 
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entropy is A = 1 (bit) because the eigenvalues of p(0) are {0, 0, 1/2, 1/2}. As time goes on, 
the entropy becomes larger because the thermalization makes the system reside equally in 
all states. Eventually, the entropy will reach to the maximum value of A = 2 (bits) where 
all four states are equally probable. 

To examine the effect of the perturbing field on an entangled state, we now consider a 
different initial density operator. We assume that the system is in a pure spin singlet at 
t = 0; \ip ) = (| t|) — | |t))/v2 and its density operator is p(0) = l^oX^ol- In Fig. 2- 
(a), we plot the diagonal components of the density operators in the multiplet basis as a 
function of time. P44 (solid line) loses its coherence linearly to time while other components 
Paa grows as time elapses. This behavior gives rise to an increasing value of the entropy 
as shown in Fig. 2-(d). For the pure initial state one can calculate the fidelity of the gate 
without too much difficulty We compare fidelity ((ipo\Ul{i)p(t)\ipo)) of a given entangled 
pure state (dotted line) with the gate fidelity (solid line) in Fig. 2-(c) and in addition the 
purity (trp(t) 2 ) of a given initial entangled state (dotted line) with gate purity (solid line) in 
Fig. 2-(d). In both quantities, there are a slight difference between the cases. This implies 
that although the gate fidelity T and gate purity V define the global characteristics of gate, 
fidelity and purity of the gate for a specific input state depends on input itself. 

Now, we discuss the strength of the decoherence which depends on T(t) and A(t) of Eq. 

\ 2 fi r^ c [a; 1 X 2 i] f UJc 

T(t)+iA(t) = / uj cos ut coth — duo — i / u>smu>tdu>. (54) 

7T jo \_2kBT\ n Jo 

For a sufficiently high temperature ksT ^> u c /2, T(t) and A(t) are further simplified to 

" sin u r t cos uj r t 



ru\ , -Am 2r sin^ c t A 
T(t) + iA[t) = 1 — 

7TT. t r s 



(55) 



UJ c t 2 t 

with Tq = X 2 rjkBTT s and Ao = X 2 t]uj c t s /7t. Since a typical value of t s is 25psec for Jo = 
and, thus uj c t s ^> 1, T(t) and A(t) are rapidly oscillating functions. This implies that the 
dominant contribution to the decoherence can be written as T(t) +iA(t) = 2To5(t)/T s in the 
limit of uj c t s ^> 1. In this approximation, we find that p a /3\-y5{t) of Eq. ( [501) is proportional to 
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Tot. This behavior is attributed to a linear dependence of various quantities(s, J 7 , V) on time. 
In addition, we expect that the degradation of the spin polarization is also proportional to 
Tot. For this, we examine the evolution of the spin polarization of the first electron for the 
initial density operator of Eq. ([51]) for various values of the coupling constant, \ 2 rj, and plot 
results in Fig. 3- (a). As X 2 r] increases, we find that more strong decoherence occurs in spin 
states and its dependence is linear on X 2 r] as shown in Fig. 3-(b). This linear dependence 
also appears in the fidelity and purity. 

In summary, we first derive an exact reduced-density-operator for the output quantum 
states in time-convolutionless form by solving the quantum Liouville equation for a noisy 
quantum channel. The formalism developed in this paper would be general enough to 
model a noisy quantum channel if various Hamiltonians for a channel dynamics, environment 
and an interaction are given. Secondly, we calculated various characteristics including the 
fidelity, purity, and the change of entropy of a two-bit quantum gate which is based on the 
spin exchange interaction between two quantum dots. Our calculation shows it is really 
important to control the decoherence in the quantum gate to protect quantum information 
against corruption. The decoherence in the quantum logic gate which is extremely sensitive 
to it may be a major obstacle to building the realistic quantum computer, however, it 
it is expected that as long as the error rate is below some threshold value, a quantum 
computer which can give arbitrary accurate answer can be built with a reasonable model 
of decoherence. In this respect, it will be interesting to investigate the implementation of 
quantum error correction technique for this model. Another interesting study on the present 
model is to find an operator sum representation for the evolution operator 8: 

£[p}=J2A,pAl (56) 

where is an operator acting on the system alone. With the operator sum representation, 
we can calculate various information theoretical quantities such as the coherent information, 
entropy exchange, and the channel capacity PT5| . We would like to leave this subject for 



future work. 
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Figure Captions 



Fig. 1. The calculated spin polarization(s), fidelity (JF), purity ("P), and entropy (A) are 
plotted as a function of time (solid lines), and compared with those obtained in the Ref. fl5( 
(dotted lines). We assume that the first electron is un-polarized on the initial state with the 
second polarized upward. For < t < t s , the swap operation is made by turning on J(t) 
and, then, J(t) = 0. 



Fig. 2. For the initial density operator p 44 , we show diagonal components of the density 
operator as time elapses in (a), pa decreases monotonically(solid line) whereas others of 
diagonal components become larger (dotted lines). In (b) and (c), the fidelity and purity are 
shown concerning with(solid line) and without (dotted line) an average over initial states. 
The evolution of the entropy as plotted in (d) starts from zero because the initial state is 
pure. 



Fig. 3. For various values of the coupling constants, X 2 r] = 0.5 x 10~ 5 (dotted), 1.8 x 
10 -5 (solid), and 3.0 x 10~ 5 (dashed), we show the evolution of the spin polarization of the 
first electron in (a) for the initial density of Eq. (|5lD . In (b), the degradation of the spin 
polarization(s), fidelity (.F), purity (V), and entropy(A) are compared for different parameters 
X 2 rj at t — t 8 . 
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Fig. 3 



